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Abstract 

The paper describes the occurrence of a pair of oppositely directed sudden impulses 
( SI), in the geomagnetic field ( AX), at ground stations, called SI + - SI~ pairs, that 
occurred between 1835 UT and 2300 UT on 23 April 1998. The SI + - SI~ pair, was 
closely correlated with corresponding variations in the solar wind density, while 
solar wind velocity and the southward component of the interplanetary magnetic 
field (Bz) did not show any correspondence. Further, this event had no source on the 
visible solar disk. However, a rear-side partial halo coronal mass ejection ( CME) 
and an Ml. 4 class solar flare behind the west limb, took place on 20 April 1998, 
the date corresponding to the traceback location of the solar wind flows. This event 
presents empirical evidence, which to our knowledge, is the best convincing evidence 
for the association of specific solar events to the observations of an SI + - SI~ pair. 
In addition, it shows that it is possible for a rear side solar flare to propagate a shock 
towards the earth. 

Introduction 



It is well known that space weather events observed at 1 AU are all linked to the 
dynamic evolution of the solar photospheric magnetic field. This evolution, in con- 
junction with solar rotation, drives space weather through the continuously changing 
conditions of the solar wind and the interplanetary magnetic field (IMF) within it. 
In spite of the fact that there have been substantial observations and discussions 
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on the close correspondence betwee n solar wind p arameters at 1 AU and groun d 
based geomagnetic fie l d variations lounsev ( 1961 ): Heppner and Mayna~ra\ (1987); 
Goodrich et all (| 19981) : \Lu et all (| 1998b it is not straightforward, under this broad 
framework, to pinpoint either the solar origins of specific space weather events or 
find specific correlations between solar wind parameters at 1 AU and ground based 
magnetic observations. This is because such signatures are generally weak and are 
usually washed out or masked by the large variety of interactions that can take place 
both in the interplanetary medium and within the earth's magnetosphere. Space 
weather events are however, often preceded by the arrival at 1 AU of strong inter- 
planetary (IP) shocks. Since such storms can have adverse effects on human tech- 
nologies, the study of IP shocks can yield important inputs for numerical models that 
simulate the propagation of solar-initiated IP disturbances out to 1 AU and beyond. 

On the other hand, solar sources of space weather events can range from coronal 
mass ejections (CME), very energetic solar flares, filament eruptions and corotating 
interaction regions (CIR). Though a vast majority of such events are caused by ex- 
plosive and energetic solar events like CME's and flares, some recent studies have 
unambiguously associated large space weather events at 1 AU, like "solar wind dis- 
appearance events", to small transien t mid-latitude corona holes but ting up against 
large active regions at central meridian \janardhan et al. ( 2005 . 20Q8al bh. These stud- 
ies have provided the first observational link between the sun and space weather 
effects at 1 AU, arising entirely from non-explosive solar events. 

Though the very first observations, by the Mariner 2 spacecraft in 1962, of in- 



terplanetary shock waves showed the po ssibility of 



ensembles in the interplanetary medium \Sonett et al. 



t he exi stence of double shock 
(| 1964b . the existence of such 



shock pairs was firmly established only some years later, by the c areful analysis o f 
plasma and magneti c field measurements associated with shocks IBurlagal (|1970b : 



Lazarus et all (|1970|) . However, the very unusual plasma and field variations asso- 



ciated with these structures prompted \Sonett and Colburm (| 1965b to suggest that the 
first or 'forward' shock would give rise to a positive H-component at ground based 
observatories while the second or 'reverse' shock would cause an oppositely directed 
or negative change in the H-component of the earth's horizontal field, as measured 
along the local geomagnetic meridian (H). These impulses, referred to in the rest of 
the paper as sudden impulse or SI + - SI~ pairs, were typically separated by a few 
hours in time and were hypothesized, as already stated, to be caused by the arrival 
at 1 AU of t he forward and reve rse shock pair convected towards the earth by the 
solar wind. \Raz.dan et all (| 1965b described worldwide occurrences of such SI + - 
SI~ pairs and suggested that they were associated with solar disturbances driving 
interplanetary shocks at highly oblique angles to the solar wind streaming direction. 
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They however, did not find any solar activity or associated occurrences of solar radio 
emission during the period of SI + - SI~ pairs. In a more recent study o f a number 
of SI + - SI~ pairs, covering the period 1995-1999. uakeuchi et all (|2002|) concluded 
that the observed SI~ (or negative impulses) in their sample were not associated with 
reverse shocks and showed no pre- ferential association to any particular kind of so- 
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Figure l : The first six panels starting from the top show the variations as a function of time in UT on 
23-24 April 1998 of the solar wind velocity, Bz, the Electric field, the solar wind flow pressure, the 
solar wind density (shaded grey) and the the symmetrical H Field respectively. Observations of the H 
field at Indian geomagnetic stations Gulmarg, Alibag and Trivandrum are shown in the bottom most 
panel. The pair of dashed vertically oriented parallel lines in all panels demarcate the times 1835 
UT and 2300 UT. These times correspond respectively to the times of the SI + impulse and the SI~ 
impulse in the H field, that were observed at Indian stations. 



lar wind structure, like high and low speed stream interface discontinuities or front 
boundaries of interplanetary magnetic clo uds. 



Early theoretical support came from \Dryen (1197(1 Il972r) who introduced the 
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physics of finite electrical conductivity within the forward and reverse shock pairs, 
to derive reasonable first order predictions for the observed distribution of solar 
wind speed , density, and temperatur e and this was fo l lowed up by severa l other 
early paper sEviatar and Dryen (|1970l ): \Shen and Dryer d 1 972b : hryer et all (|l975h 
on similar lines. In more recent times, there have been a number of theoretical 
models that have used in puts from s ol ar data to predict the ar rival of IP shock s 



and IP CMEs at eart h \odstrcd (|2003h : Wandegriff etaA (|2005b : \T6th et all d2005h : 



(HAFV2), \Fryet al. 



Detman et al\ (12006 ). including the well known Hakamada-Akasofu-Fry model 



(|2003r) which is the only model, to da te, to have been su b 



stantially validated in an operational forecasting environment 
during solar cycle 23. 



Smith et all (|2009aU bh 



The SI+ - SI Pair of 23 April 1998 

An SI + - SI~ pair was identified at three Indian geomagnetic observatories on 23-24 
April 1998. Figure Q] (bottom-most panel) shows the tracings of H magnetograms 
(projected onto the X or geographic north direction and marked AX in Fig. [[]) on 
23-24 April 1998 at the three Indian stations Gulmarg, Alibag and Trivandrum re- 
spectively. A sudden positive impulse in H was recorded at all three Indian obser- 
vatories at 1835 UT (23 Apr. 2335 LT) followed by a sudden negative impulse at 
2300 UT (24 Apr. 0430 LT). During the time interval between the SI + impulse and 
the SI~ impulse, the amplitude of H first decreased and then attained a peak of 44 
nT at Trivandrum that progressively increased to 54 nT at Alibag and 76 nT at Gul- 
marg. Between 2100-2300 UT large fluctuations were recorded at all stations. The 
fluctuation in H, at all stations were remarkably similar to each other with the am- 
plitude increasing from Trivandrum to Gulmarg. Also shown in FigH] (starting from 
the top and going down) are the corresponding variations of the solar wind velocity, 
the IMF-Bz, the interplanetary electric field, the solar wind flow pressure, the so- 
lar wind density and the symmetrical H Field respectively. The curve f or the solar 



wind density, as observed by the Advanced Composition Explorer (ACE; Stone et al. 



(| 1998b ) has been shaded grey, in FigJH in the region between the SI + - SI~ pair. The 
vertically oriented dashed parallel lines in all panels demarcates the time interval be- 
tween the SI + impulse and the SI~ impulse. It is to be noted that the symmetric H 
field (S YM/H), characterizing the mean variation of H at all middle latitude stations 
around the world, too had remarkably similar variations as the H at Indian stations. 
Thi s therefore implies that t he SI + - SI~ pair was a global event. 

Rastogi and Patel ( 1975b had shown that solar plasma moving towards the earth's 
magnetosphere with the velocity, V and having a frozen-in magnetic field normal 
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to the ecliptic (IMF-Bz) is equivalent to an electric field, E sw = (-V x Bz), which 
is transmitted without any time delay to the polar region and then to the low lati- 
tude ionosphere. This belongs to a process known as o yershield i ng ele c tric field in 
the magnetosphere which has been extensive l y studied NishidA ( 1968 ): Vasyliunas 
dl970b : \Spiro et ali (|l988h : \wolfet ali fepOlh : \Goldstein et ali ilOOj . Prompt pen- 
etration occurs due to the slow response of the shielding electric field at the inner 
edge of the ring current that opposes the time varying convection field, in the pres- 
ence of an IMF-Bz. The time scale of thi s proce s s is generally of th e order of an hour 
but can sometimes be longer I Vasyliunas (|1970h : \Senior and Blana (| 19841 ). During a 
period of sudden northward turning of the IMF, from a steady southward configu- 
ration, the convective electric field shrinks while the shielding electric field takes a 
longer time to decay and produce a residual electric field, known as the overshield- 
ing electric field. The direction of this field is opposite to the normal direction of 
the ionospheric electric field. The E sw has a direction of dusk-to-dawn for positive 
IMF-Bz and dawn-to-dusk for negative IMF-Bz. 

It can be seen from Fig. [I] that AX at Indian stations just after the SI + at 1835 
UT (around local midnight) had gradually decreased till 2000 UT. This effect is due 
to the prompt penetration of the electric field when the IMF-Bz is negative during 
the period. At around 2000 UT, AX again increased suddenly to values much above 
the first impulse level. Correspondingly, it can be seen that IMF-Bz turned from 
southward to northward at this instant implying that it is due to the overshielding 
condition described above. After this, the level of AX went down, with some oscil- 
lations and finally came down to normal level at 2300 UT. It is interesting to note 
that the fluctuations in AX between 2130-2200 UT were very well correlated with 
the solar wind density rather than with the IMF-Bz or the solar wind speed. The SI + 
at 1835 UT was associated with a sudden increase of both the solar wind densit y and 
speed causing a sudden pressure on the magnetosphere (as first suggested by Gold 
(|1959h ). The SI~ at 2300 UT was associated with the sudden decrease of solar wind 
density. 



Global Geomagnetic Fields 

Figure |2] shows the variation of AX from eleven low latitude stations around the 
world on 23-24 April 1998. The stations, starting from Alibag (ABG - Lat. 18.64; 
Long.72.87), (uppermost curve) and arranged in increasing order of geographic lon- 
gitude are respectively, Gnangara (GNA - Lat. -31.78; Long. 115.95), Esashi (ESA 
-Lat. 39.24; Long. 141.35), Honululu (HON - Lat. 21.32; Long. 202.00), Frenso 
(FRN - 37.10; Long. 240.30), Del Rio (DLR - Lat. 29.49; Long. 259.08), Kourou 
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Figure 2: Variation of the horizontal component of the geomagnetic field projected on to the x- 
direction at eleven low latitude stations around the world on 23 - 24 April 1998. The vertical dashed 
line is marked at 20.5 UT and shown alongside it, for each curve, is the corresponding local time at 
each station. Also indicated on the left of each curve is the geographic longitude of the station. 



(KOU - Lat. 2.21; Long. 307.27), Ascension Island (ASC - Lat. -7.95; Long. 
345.62), Hermanus (HER - Lat. -34.42; Long. 19.23), Addis Ababa (AAE - Lat. 
9.02; 38.77) and Tanananarive (TAN - Lat. -18.92; Long. 47.55). The geographic 
longitude of the ground stations are indicated at the left of each curve in Fig. |2l Also 
indicated to the right of the vertical dashed line (marked at 20.5 Hrs UT in Fig. 0) 
for each curve, is the local time at 20.5 hrs UT. The stations chosen range from ge- 
ographic longitudes of 19° to 346° corresponding to local times of ~22 hr. through 
the midnight, dawn, noon to dusk (20 hr). 

The negative IMF-Bz between 1835 - 2000 UT caused a decrease of AX at sta- 
tions in the night sectors (HER, AAE, and TAN ) and an increase at stations in 
the mid day sector (FRN, DLR, KOU and ASC). Around 20.50 UT, AX showed 
strong positive peaks at AAE, TAN and ABG, no change at ESA and HON and neg- 
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ative peaks at FRN, DLR and KOU. This data conform s very w e ll with the process 



Vasvliunas 



of prompt penetration and over s hielding electric fi eld Nishida f 1968h : _ 
(|1970i) : \Rastogi and Patea (| 19751 ): \Spiro et al. wherein, a southward IMF-Bz 

(between 1835 and 2000 UT) would cause a decrease of AX at night side stations and 
a increase of AX at day side stations of the earth while a northward turning of IMF- 
Bz would produce a strong positive AX at stations in the night sector and negative 
AX at stations in the day side sector due to the imposition of either a dusk-to-dawn 
or dawn-to-dusk electric field. The fluctuations in AX between 2130 - 2300 UT are 
synchronous at all stations in the day as well as in night sectors, suggesting the effect 
to be due to solar wind flow pressure and not due to the IMF-Bz. It is important to 
note here that the solar wind density fluctuations virtually mirrors those seen in AX 
at ground stations, thereby implying that the solar wind density was the main key 
or driver for this event. Qualitatively, the fluctuations seem to be independent of the 
latitude. The dominant parameter is the solar wind pressure that makes the magne- 
tosphere shrink and expand self similarly, with some scaling factor depending on the 
pressure. This is reflected in the magnetic field data at all latitudes and longitudes 
in a configuration where the IMF appears to have no role to play. Thus, this was a 
unique space weather event in which one could unambiguously associate solar wind 
density variations with variations in AX at ground stations while no such changes 
were seen in the solar wind speed or magnetic field. 



Solar Source Locations 



It is well known that due to the rotation of the sun (Q. = 1.642 x 10~ 4 deg s -1 ), a radi- 
ally directed outflow of solar wind from the sun will trace out an Archimedean spiral 
through the interplanetary medium. For a steady state solar wind with a velocity of 
430 km s _1 , the tangent to this spiral, at 1 AU, will make an angle of 45° with the ra- 
dial vector from the Sun \Schwenn\ (|l9901 ). As a consequence, the longitudinal offset 
((f) r) of a solar wind stream with a velocity v, when traced backwards from a distance 
Ri (say 1 AU) to a distance R2, will be (f)R = fl(Ri - R2)/v. We can thus project the 
observed solar wind velocities back to the sun to determine the sources of the solar 



wind flows at the sun. The earliest instance o 



using such a technique to trace solar 



wind outflows back to the sun was by \Ricked (| 19751 ). For the present event, we have 
back-projected the observed ACE velocities along Archimedean spirals to the source 
surface at 2.5 R to determine its solar source location. Though this method is gen- 
erally applicable to a steady-state flow of the solar wind, it has also been applied 
in cases when the solar wind outflows were not steady-state and highly non-radial. 
For example, during the well know disappearance event of 11 May 1999, the work 
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by Uanardhan et al\ (120051 ): Uanardham (|2006l ) has shown that solar source locations 
determined by the traceback technique, using constant velocities along Archimedean 
spirals, do not have significant errors even though the solar wind flows were known 
to be highly non-radial during that period. In the case under discussion the solar 
wind flows would have been highly kinked and non-radial due to the propagating 
forward and reverse shocks arising from the optically occulted flare and the rear side 
CME. Therefore, if the SI + - SI~ pair had a source on the solar disk, the ambiguity 
abo ut the location o f the s o urce region would be within reasonable errors as shown 
by Uanardhan et al.\ (|2005h : Uanardham (|2006l ). Figure [3] shows a map of the solar 
photosphere indicating the locations of the active regions. The back projected region 
of the solar wind flows go back to the vicinity of the large active region AR8205 



20 April 1998 




Figure 3: Map of the solar photosphere on 20 April 1998 corresponding to the back projected region 
of the solar wind flows obsrved at 1 AU. The map shows the locations of the large active regions with 
active region AR8205 indicated by an arrow for convenience. 

located at N21W25, to the west of central meridian on 20 April 1998, and indicated 
by a solid arrow in Fig. [H Typically, the solar disk shows a large number of ac- 
tive regions during the rising phase of the solar cycle. A detailed theoretical study by 



Schrijver and DeRosd (|2003|),has shown that solar wind outflows from active regions 
comprise <10% during solar minimum and up to 30-50% during solar maximum. 
However, the visible solar disk on 20 April 1998 showed no activity in terms of flares 
or CME's and had only two active regions AR8206 and AR8205 as shown in Figure 
[H The Active region AR8206 was smaller than AR8205 being around 245 millionth 
of the solar disk in size as compared to 312 millionth of the solar disk for AR8205. 
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Also, AR8205 was less than 30° west of central meridian as compared to over 40° 
east of central meridian for AR8206. It may be noted that a central meridian location 
would imply that any activity like a large CME or flare would be earth directed. How- 
ever, there was no flare or CME on the entire visible s olar disk on 20 April. Images 
from the Extreme-ultraviolet Imaging Teles cope (EIT: \Delaboudiniere et all dl995|)) 
and the Michelson Doppler Imager (M DI. \Scherrer et all (11995b ) onboard the Solar 
and Heliospheric Observatory (SoHO iDomingo et all (| 1995b ) were also examined 
carefully to confirm that there was no other possible source regions on the solar disk 
on 20 April 1998. 



The Rear side CME and Optically Occulted Flare of 20 April 1998 

On 20 April 1998, a rear side, fast (~1850 km s -1 ) partial halo CME occurred in 
association with an optically occulted GOES Ml. 4 class flare which took place just 
behind the limb at S43W90. It must be pointed out here that most forecasters of 
space weather events generally ignore the possibility that a limb or backside solar 
explosive event could propagate a disturbance towards the earth. However, there 

p een stu died and reported in recent 
' (|2009bfl) . 



have been some instances where such cases have 
times 



Mc Kenna-Lawlor et al\ (|2006l) ; \Smith et al. 



The GOES M1.4 flare at S43W90 was first detected in 1-8 A band at 09:15 UT 
on 20 April 1998 and reached its maximum at 10:21 UT. The rear side partial halo 
CME was first seen in the LASCO coronograph C2, at 10:04:51 UT on 20 April 
1998, as a bright, sharply defined loop structure spanning ~80° in latitude and ex- 
tending to ~3.1 R.0. The same was first observed by C3 at 1 0:45:22 UT. Both th e 
CME and the flare h ave been extensively studied and reported \Bastian et all (|200l|) : 



Gopalswamy 



Simnett (2000, 2002) and it has been shown that the CME, which was radio loud 



(|2000b . actually pushed aside pre-existing streamers while moving be- 



yond the LASCO C3 field of view. Since this was a rear-side CME the shock front 
that it drove would have been in a direction away from t he earth. In a s t udy o f 
the arrival time of flare driven shocks at 1 AU and beyond \Smart and Shea (II 985b ■ 



Janardhan et all (LL996) it has been assumed that the trailing edges of flare driven 



shock waves travel at roughly half the velocity of the shock in the flare radial di- 
rection. It is therefore not unreasonable to assume that the trailing edges of the 
CME driven reverse shock would be much slower and could be convected outwards 
towards the earth by the solar wind. The flare and the rear-side CME would thus pro- 
vide the forward and reverse shocks to cause the SI + and SI~ pair. Figure @] shows 
hourly averaged value of the absolute magnitude of the magnetic field, as observed 
by the ACE spacecraft, as a function of time in UT. It is expected that the strength of 
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the magnetic field would increase at the forward shock or SI + impulse and decrease 
at the reverse shock or SI~ impulse. It can be easily seen from Fig. |4] that there 
is a sharp increase in the magnetic field at around 1835 UT, corresponding to the 
arrival of the forward shock associated with the SI + impulse and a decrease in the 
magnetic field at 2300 UT corresponding to the arrival of the reverse shock associ- 
ated with the SI~ impulse. The vertically oriented dashed parallel lines in Fig. @]are 
marked at 1835 UT and 2300 UT, the time corresponding to the SI + and SI~ impulse 
respectively. 

20 
15 
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Figure 4: Hourly averaged total magnetic field as a function of time in UT as observed by the ACE 
spacecraft located at the LI Lagrangian point at 1 AU. The vertically oriented dashed parallel lines at 
1835 UT and 2300 UT correspond to the time of the SI + and SI impulse respectively. 

Discussion and Conclusions 

From an observational point of view, the present work has been able to link interplan- 
etary structure during this particular event with worldwide magnetospheric response, 
using the Indian magnetic observatories to provide the first clue. In particular, this 
event has been unique in that the solar wind density variations have played a key 
role, as seen through the close correspondence between the fluctuations in the solar 
wind densities and the AX at ground stations while no such changes were seen in the 
solar wind speed or magnetic field. Though there has been a large body of work, over 
the past four decades, that have addressed the issues concerned with forward/reverse 
shock pairs, their manifestation at 1 AU and their relation to specific solar events, we 
believe that this paper presents empirical evidence, which to our knowledge, is the 
best convincing evidence for the association of specific solar events to the observa- 
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tions of an SI + -SI~ pair. In addition, it shows that it is possible for a rear side solar 
event to propagate a shock towards the earth. 

We have seen that the SI + impulse 1835 UT was associated with a sudden in- 
crease of both the solar wind density and speed causing a sudden pressure on the 
magnetosphere while the SI~ at 2300 UT was associated with the sudden decrease 
of solar wind density. The southward IMF-Bz between 1835-2000 UT caused a de- 
crease of AX at night side stations and a increase of AX at day side stations of the 
earth due to the imposition of a sudden electric field caused by the prompt pene- 
tration of electric field to low latitudes. As stated earlier, a northward turning of 
IMF-Bz produces a strong positive AX at stations in a night sector and negative AX 
at stations in the day sight sector due to the effect of overshield electric field which is 
in a direction opposite to the normal electric field in the ionosphere. Between 2015- 
2300 UT the fluctuations in AX were similar at all stations in the day or night sectors 
and were well correlated with the fluctuations in solar wind flow pressure, reflecting 
the shrinking and expansion of the magnetopause as a result of strong solar wind 
pressure variation. 

The solar event lasting only for some 4 to 5 hours showed signatures of all mech- 
anisms involving solar - magnetosphere - ionosphere relationships. The arrival at 
1 AU of the forward and reverse shock pair associated with the SI + and SI~ re- 
spectively is clearly seen in the behaviour of the hourly averaged values of the total 
magnetic field, which shows a sharp increase at ~1835UT and a decrease at ~2300 
UT. The effect of sudden changes in the solar flow pressure due to a change of only 
the solar wind density have been clearly identified. The effect of the slowly vary- 
ing IMF-Bz has been shown to impose dusk-to-dawn or dawn-to-dusk electric field 
globally, depending on the southward or northward turning of the IMF-Bz. Though 
theoretical first order predictions for the observed distribution of solar wind speed, 
density, and temperature (as in Fig.l) durin g the propaga t ion of forward and reverse 
shock pairs were derived four decades ago Dryer ( 1970 . 1972b . the analysis of the 
event has been rewarding due to the relative quiet solar conditions prevailing as it 
allowed us to identify specific solar sources as the possible drivers of the SI + and 
SI~ pair. The only activity on Sun was the rear side CME and the associated solar 
flare. This is thus a very unique observation wherein, a pair of SI events have been 
shown to be associated with corresponding changes in the solar wind density while 
no such changes are seen in the solar wind speed or magnetic field. Many more such 
events need to be observed, retrieved and studied, both from archival records and 
future observations, before a clearer understanding of the exact nature and physics 
behind suc h events is obtained. High resolution, high dynamic range radio imaging 
techniques Mercier et all ( 2006 ) can also provide useful information in this regard. 
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